function [ X ] = PspatialMap( h , Lxt , Lxe , Lye , Lyt , NIt , NJt , NJe , NIe , dimP)
%SPATIALMAP Summary of this function goes here
%   Detailed explanation goes here

x = zeros (dimP,1);
y = zeros (dimP,1);

xStart = h/2;
yStart = h/2;

for i=1:NIt
    for j=1:NJt
        
        if (i<=NIe && j>NJe) continue; end
        
        index = Pmap (i,j);
        
        x (index) = xStart + (i-1) * h; 
        y (index) = yStart + (j-1) * h;
        
    end
end

X = [x,y];

end

